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ABSTRACT 

Half-lives  of  thermal  neutron  activated  samples  were 
determined  by  direct  decay  measurements „   Transport  of 
samples  between  activation  and  counting  areas  was  accomplished 
by  employment  of  a  rapid  pneumatic  transfer  system „  Use  of  a 
single  channel  analyser  with  a  fast  chart  recorder  to  record 
scaled  counts  vs.  time,  and  analysis  of  these  data  in  a 
computer  program  are  discussed .   Results  for  nuclides 
investigated: 

Hf17&n  tj  m  (3o92  *  005)  sec0 
Hfl79m  ti  s  (l£049  ♦  004)  sec0 
Yb177m      tj  s  (6o49  t  o03)  secQ 

Our  method  has  several  advantages  over  multichannel 
devices  and  appears  to  be  capable  of  measuring  half-lives  on 
the  order  of  „5  seconds  with  approximately  1%   precision 

The  authors  appreciate  the  assistance  rendered  by  the 
following  professors  in  various  phases  of  the  programs 
E.  A0  Milne  and  If.  W.  Hawes,  for  supervising  reactor  operations 
whenever  called;   J.  R*  Borstings,  for  advice  on  statistical 
matters;  and  H.  E.  Handler,  for  general  counseling  and 
assistance  throughout  the  program0  The  authors  also  wish  to 
acknowledge  the  assistance  Profo  G,  W0  Rodeback,  thesis  advisor? 
gave  during  this  projecto 
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I„   Introduction 

Comprehensive  investigations  of  isomeric  states  of 

isotopes  of  Hafnium  and  Ytterbium  have  been  reported  several 

times  since  the  discovery  of  nuclear  isomers  /l~5/o  The  results 

of  the  half -life  measurements  reported  by  EQ  Co  Campbell  and 

P.  F.  Fettweis  in  1959  A/  are  summarized  in  the  following 

table,  along  with  the  corresponding  reactions  used  and  the 

gamma  energies  detected . 

reaction  half -life  gamma  energy 

Ybl76  (n>tf)  Ybl77m    T£  s  6,5  sec        o10  Mev 
Hfl77  (n,tf)  Hfl7&n    Tj  s  3c5  sec        o^$    u44  Mev 
Ta17^  k  captureHf17am  Tj  _  ^   sec  g  Mev 

Hf!7^  (n,tf)  Hf179m    T^  =  lg6  sec       o22  Mev 
Of  interest  is  the  apparent  difference  in  their  results  for 
the  half -life  of  Hf^&n,  although  this  may  have  been  within 
the  limits  of  the  experimental  errors,  which  were  not  reportedo 
More  recent  work  (1962)  by  K0  F0  Alexander  and  H0  FG  Brinckmann 
/5.  indicates  half -lives  of  (4*3  -  0el)  sec  and  (lg96  t   0oS)  sec 
for  the  Hfl77,17d  (h)JT)  Hf 173m, 179m  tactions. 

The  objective  of  this  work  has  been  to  improve  the  precision 
of  the  half-life  measurements  of  Yb177m>  Hf178m,  and  Hfl79m0 
Measurements  were  made  by  scintillation  crystal  detection  of 
direct  gamma  decay  of  thermal  neutron  activated  samples 0  The 
time  delay  between  activation  and  measurement  was  minimized 
through  the  use  of  a  rapid  pneumatic  transfer  system  which 
transfers  the  sample  from  the  activation  region  to  the  counting 
area,  Use  of  a  single  channel  analyser  and  fast  scaler  kept 


dead  times  below  those  commonly  encountered  with  multichannel 
instruments •   A  fast  continuous-chart  recorder  was  employed 
for  accurate  recording  of  data0 

Data  was  analysed  in  a  computer  program  (FRANTIC  II ^  by 
P.  Co  Rogers  of  M0I.T.)  /6/  which  performs  an  iterative  least 
squares  analysis  of  multiple  exponential  decays «,   Several  runs 
have  been  made  for  each  nuclide,  not  only  to  check  reproducibility 
of  the  measurements,  but  also  to  determine  mean  half -life 
values  based  on  statistical  weighting  and  the  corresponding 
error  estimates  for  these  mean  half -lives • 


II 0   Experimental  Equipment  (Schematic ,  Fig0  1) 
A„   Transfer  System 

lo   The  rapid  pneumatic  transfer  system  was  designed 
and  built  for  use  in  the  AGN-201*  reactor  at  USNPGS  by 
Paas  and  Sullivan  ff/.      It  transports  a  sample,  which  is 
placed  in  a  small  light  container ,  from  the  counting  area8 
about  eight  feet  from  the  reactor  shielding ,  into  a  high 
neutron  flux  region  in  the  reactor  core0   Here  the  sample 
is  activated  by  slow  neutrons  and  then  propelled  under 
positive  pressure  back  into  the  counting  position 0   Transit 
time  from  core  to  counting  position  is  approximately  one- 
tenth  second.   The  transfer  system  is  completely  closed 9 
allowing  recycling  of  the  sample  without  disturbing  the 
detection  arrangement • 

20   The  sample  container ,  or  "rabbit",  is  made  of 
balsa  wood  in  the  shape  of  a  bulleto   It  is  4o5  cm  long, 
has  a  1„9  cm  diameter,  and  weighs  about  Id  g  including 
the  sample .   The  isotopically  enriched  sample 9    50  mg  in 
powdered  form,  was  placed  in  a  gelatin  medicinal  capsule 
which  was  then  secured  within  the  "rabbit" „   Irradiation 
and  counting  of  a  dummy  "rabbit"  (less  sample)  indicated 
no  increase  over  background  rates 0 

*The  AGN-201  is  a  graphite  and  polyethylene- 
moderated  research  reactor  capable  *?   producing  «  neutron 
flux  up  to  4  x  10-1-0  n/cm2sec  at  the  AEC  license  limit  of 
1000  watts . 
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Figure  \ ■ 
Counting  System  and  Pneumatic  Transfer  System 


3.  Samples  were  procured  from  Oak  Ridge  National 
Laboratory  and  are  highly  enriched  in  the  nuclides  under 
investigation.   Isotopic  analyses  are  included  as 
Appendix  III. 
B.   Detector  System  (Fig.  2) 

1.  A  2n  x  2"   Nal(Tl)  crystal  was  used  as  the 
scintillation  detector.  The  crystal  used  was  a  Harshaw 
Model  8S8/2-I  integral  assembly  incorporating  a  Dumont 
6292  photomultiplier  tube  with  an  S-ll  response 0   The 
assembly  is  rated  by  the  manufacturer  to  be  capable  of 
7.6$  resolution.*  Our  measurements  showed  the  resolution 
to  be  6,0$  at  the  recommended  1000  volts  anode=cathode 
operating  voltage.   During  counting,  the  crystal  was  seated 
flush  against  the  rabbit  tube  in  the  shielded  counting 
area.  The  shielding  material  was  lead  with  a  minimum 
thickness  of  two  inches.  The  counting  area  was  off -set 
about  16  inches  from  the  rabbit  tube  center-line  to  reduce 
the  background  rate  at  high  reactor  powers 0   A  typical 
value  of  background  is  one  count  per  second  at  an  operating 
power  of  20  watts. 

2.  The  crystal-PM  tube  assembly  was  connected  to 
a  6AK5  cathode-follower  preamplifier  which  relayed  the 


The  resolution  is  defined  as  the  ratio  of  the  energy 
width  at  half  maximum  to  the  mean  energy  of  the  0662  Mev 
C3137  photopeak. 


si<mal  through  50  feet  of  shielded  cable  to  the  input 
of  the  linear  amplifier 

3 »   High  voltage  for  the  photomultiplier  tube  was 
furnished  by  a  regulated  power  supply  with  a  specified 
stability  of  less  than  .0026$  output  change  per  line  volt 
input  change,  and  a  maximum  output  change  of  .2$  per  day0 
No  significant  H.  V.  drift  was  observed  during  the  decay 
runs.   Preamplifier  power  (t300v)  was  supplied  by  a 
stabilized  power  supply  unit  in  the  single  channel 
analyser. 

C.  Single  Channel  Analyser 

A  commercial  combination  nonoverload  linear  amplifier 
and  integral-differential  discriminator  was  usedo 
Continuous  amplification  (voltage  gain)  from  unity  to  960 
was  available  with  the  linear  amplifier.  Amplifier 
linearity  (section  III-A)  and  an  acceptable  signal-to-noise 
ratio  (45:1  for  Cs-^-37  photopeak)  were  verified  for  gain 
settings  in  the  range  used  for  all  runs0  The  discriminator 
unit  allowed  continuous  selection  of  pulse  heights  from 
0-100  volts  (normal  calibration  units)  and  differential 
pulse  selection  from  0-10  volts.  The  whole  assembly  had  a 
resolving  time  of  approximately  one  microsecond. 

D.  Count  Recording  System  (Fig.  3) 

1.  The  output  of  the  discriminator ,  a  20  volt  pulse 8 
was  transmitted  to  a  decade  scaler  having  better  than  one 


s 

1 

I 

1 

microsecond  resolving  time0  The  scaler  has  been  modified 
by  the  addition  of  a  6CB7  cathode  follower  which  acts  as 
an  impedance  matching  device „  Thus  when  the  scaler  was 
operating  in  a  continuous  mode,  scaled  output  pulses  could 
be  fed  to  the  chart  recorder  without  stopping  or  recycling 
the  scaler.   Scaled  output  pulses  could  be  selected  for 
each  2,  10,  102,  or  103  input  pulses  to  the  scaler 0  The 
scaled  pulses  were  sent  through  a  differentiating- clipping 
circuit  and  then  to  the  continuous  chart  recorder0 

2.  The  continuous  recorder  is  a  two-channel  chart 
type  with  variable  chart  speeds  of  ls  5»   25,  and  125  mm/sec 0 
It  is  capable  of  recording  up  to  about  110  cps  with  a 
minimum  input  sensitivity  of  10  m0vo   For  this  application 9 
output  pulses  from  the  scaler,  representing  a  scaled  number 
of  counts,  were  recorded  in  one  channel  and  a  calibrating 
time  mark  was  recorded  in  the  other0 
E.   Time  Calibration  System 

The  basic  instrument  for  this  system  is  a  crystal 
controlled  oscillator  designed  primarily  as  a  frequency 
standard  for  radio  frequency  calibration 0   Frequency 
standards  of  10,  20,  or  100  kilocycles ,  with  rated  stability 
of  „005$,  could  be  selected 0   Pulses  of  desired  frequency 
were  transmitted  via  shielded  cable  to  a  five-decade,  200 
kilocycle  scaler.  The  scaler  was  modified  to  act  as  a 
frequency  divider  by  tapping  off  the  third  or  fourth  stages, 


thus  giving  a  wide  selection  of  time  calibrating  pulses 
(10-1,  2.9   10,  etc,  seconds  per  pulse  interval )  „  The  time 
pulses,  after  differentiation  and  clipping,  gave  a  single 
sharp  spike  per  scaled  number  of  input  cycles  and  an 
accurate  record  of  time  on  the  recorder  chart „ 
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III.  Calibration  and  Procedure 
A.   Energy  Calibration 

1.  A  planchette  containing  an  evaporated  Cs  -*'  sample 
was  used  as  a  calibrating  source .  The  known  ^  energy 
(662  Kev)  and  long  half  life  (30  yr. )  suited  it  to  this 
use.  The  planchette  was  secured  in  a  phenolic  holder  which 
was  machined  to  fit  over  the  crystal  case  (Fig0  4)»   The 
holder  and  crystal  case  were  scribed  so  that  a  reproducible 
position  was  available. 

2.  The  amplifier  linearity  was  adjusted  according 
to  the  manufacturer^  instructions.  An  output  pulse  from 
a  mercury  switch  pulser  that  matched  the  pulse  shape  of 
the  signal  from  the  preamplifier  was  used.   At  the  same 
time,  for  convenience,  the  amplifier  was  calibrated  so 
that  the  lower  discriminator  (E)  settings  gave  the 
approximate  pulse  height  in  volts.  The  use  of  Cs1^'  for 
an  energy  calibration  source  and  knowledge  that  the 
primary  gamma  energies  of  the  nuclides  of  interest  lay 
under  662  Kev,  led  to  the  choice  of  approximately  90  v* 

(.9  of  full  scale)  for  the  position  of  the  Cs  ''   photoline0 
Gain  settings  on  the  amplifier  were  adjusted  so  that  the 
photoline  fell  at  roughly  90  v.  when  observed  on  a  cathode 
ray  oscilliscope.  Exact  position  of  the  photopeak  was 
determined  from  a  differential  energy  spectrum  and  found 
to  be  at  9#.5  volts  as  read  on  the  lower  level  discriminator 
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scale.  Thus  the  discriminator  dial  setting  was  calibrated 
versus  gamma  energy.  A  check  on  the  linearity  of  the 
entire  system  and  verification  of  the  calibration 
described  above  was  obtained  by  measuring  the  photopeak 
position  of  the  364  kev  gamma  from  1^31 ,  Results  are 


shown  in  Fig,  5. 


■20 iro 60   80"  100 

Discrimator  Setting  (e)  in  volts 
Fig.  5.  Energy  Calibration  Curve 

Periodic  checks  using  the  mercury  switch  pulser  as  the 
input,  and  an  oscilliscope  to  determine  pulse  heights  of 
the  amplifier  output,  indicated  no  departure  from  amplifier 
linearity  within  the  uncertainties  of  the  measurements 0 

Calibration  of  the  differential  dis criminatory  or 
"window",  gave  a  maximum  window  width  of  #,3  volts , 
corresponding  to  approximately  a  60  kev  energy  interval 
on  the  curve  shown  in  Fig.  5» 

3.  Daily  gain  adjustments  to  provide  energy  calibration 
were  made  by  a  "quick"  method  suggested  by  H0  E0  Handler 
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Initially,  this  method  required  the  selection  of  a 
reference  discrimination  setting  and  measurement  of  the 
corresponding  integral  rate  immediately  after  obtaining 
the  Csl37  spectrum,  A  value  near  the  midpoint  of  the 


dn 
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reference  setting. 
Compton  Edge 


(36,500  t   190)cpm 


Vbv  '89.5v 


—  dE 


Discriminator  (E)  Setting 
Fig.  6.   Csl3?  Spectrum  Showing  Reference 
Voltage  and  Reference  Rate 

"valley"  between  the  Compton  edge  and  the  photopeak  was 
chosen  as  the  reference  discriminator  settingo  This 
choice  of  setting  makes  the  corresponding  integral  count 
least  susceptible  to  small  gain  changes .   At  this  setting^ 
76  volts  on  the  lower  discriminator ,  the  integral  rate 
(corresponding  to  the  area  under  the  curve)  was  measured 
as  (36,500  t   190)  counts  per  minute .  Once  this  rate  was 
determined  it  became  a  reference  for  future  measurements 0 
Subsequent  energy  calibrations  were  made  by  first  adjusting 
the  amplifier  gain  to  give  the  approximate  desired  pulse 
height  of  the  reference  photopeak  as  measured  on  an 
oscilliscope.  Then  fine  gain  adjustments  were  made  at  the 
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reference  discriminator  setting  to  attain  the  reference 

counting  rate. 

B.   Gain  Stability 

The  fact  that  certain  photomultiplier  (PM)  tubes 
exhibit  marked  gain  shifts  (effectively  a  change  in 
multiplication)  at  high  counting  rates 9   particularly  rates 
greater  than  5000  counts  per  sec  ,  has  been  reported  by 
several  investigators  /9,  10 ,  11/ .  D0  F0  Covell  and 
B.  A.  Euler  at  NRDL,  in  an  extensive  study  of  various 
brands  and  types  of  tubes,  found  that  of  the  tubes  showing 
gain  shifts,  all  required  an  hour  to  several  hours  to 
recover^  R.  Chery  corroborated  this  information  when  he 
found  PM  tubes  would  be  very  stable  for  periods  of  about 
24  hours  if  they  were  first  exposed  for  20-30  hours  to 
high  (just  below  tube  dissipation  limits)  gamma  radiation 
levels o 

This  slow  tube  recovery  has  been  utilized  to  advantage 
in  this  experiment  by  presaturating  the  PM  tube  prior  to 
the  decay  runs.  The  sample  was  irradiated  at  the  same 
power  level  and  for  the  same  irradiation  time  as  for  the 
actual  decay  run.  Reference  rates  were  checked  after  this 
exposure  and  readjusted,  as  necessary0   After  decay  runs? 
reference  rate  checks  showed  less  than  105%  gain  shift 
(in  our  case  less  than  075%   reference  rate  shift),  and 

*The  initial  gain  shift  occurred  almost  instantaneously « 
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tube  recovery  times  lasted  up  to  several  hours „   In  all 
cases,  recovery  times  were  very  much  greater  than  the 
maximum  counting  period  used  (five  minutes )e  The  effect- 
iveness  of  the  technique  was  substantiated  in  subsequent 
analyses  of  the  data  (Section  IV),  which  indicated  no 
appreciable  systematic  errors  as  would  have  been  the  case 
with  large  gain  shifts 0 

An  alternate  method  for  checking  gain  shift  consisted 
of  placing  a  Cs-^37  sample  near  the  crystal  during 
measurement  runs  to  give  a  steady  high  background  rate0 
This  rate  (^"150  cps)  was  measured  before  and  after  each 
run,  at  operating  powers ,  This  method  confirmed  that  gain 
shift  was  generally  less  than  1.0$,  and  data  were  not  used 
if  the  gain  shift  exceeded  105$=  Two  advantages  were 
gained  with  this  check.  The  first  was  that  the  measure- 
ments could  be  taken  without  disturbing  the  crystal 
assembly.   Second,  this  allowed  a  visual  check  during  the 
run  ty  superimposing  the  Cs^3'  pulse  over  the  measured 
pulses  on  the  oscilloscope .   In  this  manner  any  large 
short-term  shifts  during  measurement  could  be  observed , 
although  none  were  observed 0 

Other  electronic  equipment  in  the  counting  system^ 
once  stabilized  by  a  three  to  four  hour  warm-up  period , 
had  negligible  drift  over  the  periods  of  time  required  to 
make  one  or  several  runs. 


16 


Co   Measurement  of  Resolving  Time 

A  double  pulse  generator  was  used  to  determine  the 
resolving  time  of  the  counting  system0   Initially ,  each 
component  (amplifier,  discriminator ,  and  scaler)  was 
checked  separately  to  assure  specifications  were  meto 

A  test  signal  circuit  is  provided  in  the  preamplifier 
system  that  allowed  checking  the  system  resolving  time 
with  the  double  pulse  signal 0  The  results  of  measurements 
on  the  entire  system  gave  a  value  of  (lo4  Sol)  micro-seconds 
for  the  resolving  time9 
Do   Activation  Parameters 

Activation  times  and  reactor  powers  were  varied 
according  to  the  half  life  measured  and  the  gamma  energy 
investigated o  An  effort  was  made  in  each  case  to  have 
an  initial  recorded  count  rate  of  5»000-109000  counts  per 
second o  This  limitation  was  imposed  by  the  use  of  the  100 
counts  per  mark  scaling  factor  at  the  beginning  of  each 
runc  Rates  greater  than  10 9 000  counts/sec  approached  the 
recorder  frequency  limit  of  about  110  cycles/sec 0 

The  Hafnium  activities  were  measured  at  gamma  energies 
of  (330  t   30*)  kev  for  Hf17gm  and  (215  t   30)  kev  for 
Hf179mo   BfL7&a  also  has  a  very  strong  (41%)  line  at 
213 o4  kev  /ll/,  but  the  low  activation  cross-section  ratio 


One-half  the  energy  "window",  rather  than  an  error 
estimate o 
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hi   of  O^actdyS^actd??)  s  o054  for  our  sample 
precluded  the  use  of  this  energy  for  the  measurement 
of  the  shorter-lived  nuclide 0  To  eliminate  possible 
degradation  of  the  results  for  the  longer  half-life  due 
to  pulse  pile-up  ,  the  shorter-lived  nuclide  was  allowed 
to  decay  before  measuring  the  longer. 

In  general ,  activation  times  of  approximately  one 
half-life  were  selected  in  order  to  obtain  at  least  50% 
of  the  saturation  activity  of  the  nuclide  of  interests 
yet  reduce  contributions  to  the  counting  rate  from  any 
long-lived  impurities  present „  This  method  proved  inadequate 
for  the  investigation  of  the  Hfl?£m  nuclide 0   Consequently 9 
a  higher  reactor  power  for  a  shorter  irradiation  period 
was  chosen  to  give  a  ratio  of  initial  relative  activities 
of  at  least  3:1  for  the  Hf^&n  s  HfWm  nuclides . 

The  integral  mode  of  operation  was  used  in  the 
Yb  ' 'm  measurement o  The  low  cross-section ,  about  7  barns 
for  Yb  '  ,  would  have  required  an  excessive  amount  of 
reactor  power  in  order  to  obtain  desired  counting  rates  in 
the  differential  mode0  The  integral  discriminator  (E  dial) 
was  set  to  eliminate  all  gamma  energies  below  about  95  kev0 
In  this  manner ,  both  the  104  and  212  kev  gammas  were 
measured  at  a  considerable  savings  in  reactor  power 0 

'Tin  the  case  in  which  two  half  -lives ,  say  tj"  and  tj*» 
are  represented  in  the  source,  the  effective  pulse  intro- 
duces three  false  half-lives  into  the  measured  activity* 

IS 


The  following  table  summarizes  the  activation 
parameters  and  the  gamma  energies  for  each  type  run0 
Nuclide   Activation  Time   Reactor  Power  Gamma  Energy 
Hf17Sm        2  sec        20  watts      (330  +  30)  kev 
Hfl79m       20  1  (215  *  30) 

Yd1??1*        5  10  all  above  90 


19 


IVo  Data 

A0   Data  Reduction 

Data  tapes  consisted  of  two  channels  of  simultaneously 
recorded  information;  one  of  time,  the  other  of  scaled 
counting  marks  (Fig0  7)°  The  "counts"  channel  was  broken 
down  into  a  series  of  consecutive  time  intervals,,  starting 
with  intervals  that  were  short  compared  to  the  measured 
half =life  and  increasing  in  length  (time)  as  the  counting 
rate  decreased  and  approached  the  background  level 0  Short 
intervals  were  chosen  so  that  the  sample  would  not  decay- 
appreciably  over  any  given  counting  intervale  The 
approximate  expression  for  the  measured  rates 

ft   +  £t 
Rn  s  -h  )  t  _<£   "  dt» 

is  valid  only  for  the  case  of  At  small  compared  to  t^ 
where  n  is  the  instantaneous  counting  rate  at  the  midpoint 
(t)  of  the  n  th  interval  (At)  and  Rn  is  the  measured  number 
of  counts  per  interval  At „   In  each  of  our  decay  analyses , 
intervals  of  less  than  1/7  of  the  (shortest)  half -life 
were  used  initially.,  To  ensure  adequate  statistical 
accuracy  for  all  intervals  containing  more  than  2500  counts s 
25  marks  at  a  scaling  factor  of  100 ,  the  number  of  counts 
was"  estimated  to  the  nearest  50  counts 0  Similarly ,  the 
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number,  for  intervals  containing  625^2500  counts ,  was 
estimated  to  the  nearest  25  counts „   In  this  manner  the 
counts  in  each  interval  were  estimated  to  within  the 
standard  deviation  calculated  according  to  Poisson 
statistics o 

The  timing  error  was  estimated  by  measuring  the  drift 
of  the  calibrating  marks  in  the  "time"  channel  with  respect 
to  the  printed  grido  Typical  of  the  values  obtained  in 
these  measurements  (taken  over  the  entire  counting  period) 
was  1  x  10~3  sec/sec „ 
B0   Analysis  (least  squares  program) 

5 

1 o   Theory 

The  data  for  the  two  component  decay  of  the  Hf 

sample  and  for  the  single  decay  of  the  Yb  sample  were 

analysed  with  the  aid  of  a  CDC  1604  computer0   In  the 

iterative  least  squares  analysis  program  (FRANTIC)* 

the  input  data  points  (R^)  are  corrected  for  dead 

time  (T, )  and  background  (B)Q  Then  each  point  is 
a 

normalized  (XNORM)  as  desired  (in  this  case  XNORM 
-  loOO  was  used)  to  yield  the  correct  data  (An)  for 
the  n^  data  point, 

(1)  t     .  \~-M~         -  b]  XNORM 

A  statistical  weight  factor  (Wn)  ±s   then  applied , 

This  discussion  will  necessarily  be  limited  to  the 
main  points  of  the  analysis „  For  further  clarification 
refer  to  item  #6,  Bibliography » 
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which  is  of  the  forms 
(2)  Wn 


#..2 (total )n  (XN0RM)2 

where  CX*(total)n9  the  total  variance  for  the  n**1  pointy 
consists  of  the  sum  of  the  variances  for  1)  count  rate 
uncertainty  2)  background  uncertainty  3)  uncertainty  in 
dead  time  and  4)  uncertainty  in  the  time  length  of  the 
counting  interval 0 

For  finite  observation  intervals  (DTn)s  the  expression 
for  the  activity  (AC  )  resulting  from  N  radioactive  decays 
is; 

(3)  ACn-2  Ac.  e  ^  *  1-e      _ 

lsl   l  XLDTn 

This  exact  expression  is  not  properly  evaluated  by 
the  computer  for  small  values  of  /\DTn  so  that  equation  (3) 
is  expressed  in  terms  of  a  hyperbolic  sine  function  and 
expanded  to  give  the  expression; 


N       -yTn  ♦  DTn) 
U)   ACn  -X      Ao;  e  l 

L=  1    l 


where  x 


3§  51     7? 


O  0 


In  the  least  squares  fit  analysis  this  expression  is 
used  for  values  of  X»  DTn  up  to  lo0  and  the  original 
expression  (3)  is  used  for  values  of  ^DTn  greater  than 
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loOo  The  Xj  values  are  then  held  fixed  and  the  above 
expressions  are  linear  in  the  A0  coefficients 0 

In  the  iterative  least  squares  analysis  the  exact 
expression  (3)  is  expanded  in  a  Taylor  series  about  the 
point  defined  by  the  values  of  A0°  and  X'  from  the  preceding 
iteration o 

The  iterative  process  is  continued  until  all  parameters , 
as  well  as  the  weighted  variance  of  fit  (the  square  of  the 
standard  deviation  of  the  distribution  of  residuals  about 
zero),  deviate  from  their  value  in  the  preeeding  iteration 
by  not  more  than  one  part  in  10" 8  or  until  25  iterations 
have  been  completed „ 

A  sample  of  the  computer  output  from  a  measurement  of 
Hf179m  is  included  as  Appendix  I0 

2«   Goodness  of  Fit 

Of  particular  interest  in  this  program  is  the  goodness 
of  fit  information  included  in  the  printout.  DC^8  the 
weighted  variance  of  fit  (VAR),  and  a  histogram  of  the 
distribution  of  the  residuals  (of  the  experimental  curve 
about  the  calculated  curve)  are  givenQ   In  using  the  usual 
}£  2  vSo  degrees  of  freedom  tables ,  the  weighted  variance 
of  fit  (VAR)  times  the  degrees  of  freedom  (DP)*  will  give 
a  more  realistic  measure  than  X   of  the  probability  of 


*In  this  case  degrees  of  freedom  is  defined  as  the 
number  of  input  data  points  minus  the  number  of  parameters 
allowed  to  float « 
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performing  a  better  experiment  or  finding  a  better  fit 
with  different  parameters  or  a  different  hypothesise  This 
is  true  since  the  uncertainties  in  dead  times  background 9 
and  length  of  time  interval  are  considered  in  the 
calculation  of  VAR,  whereas  they  are  not  in  calculating  X  ° 

As  an  example,  using  the  output  data  shown  in  Appendix 
I,  VAR  xDF  s  27«55S  (with  29  degrees  of  freedom) c  This 
gives  a  91 o 5$  "assurance"  that  the  calculated  fit  is  unique 
for  the  given  data,  i0e0,  that  the  hypothesis  proposed  is 
correct o  The  histogram  of  the  distribution  of  the  residuals 
is  approximately  Gaussian,  as  it  should  be0   In  this  case8 
had  a  two  component  decay  mode  been  proposed  in  the  initial 
estimates,  the  anlysis  would  have  either  rejected  the 
second  mode  or  given  a  large  VAR  indicating  a  false 
hypothesis. 

In  general,  an  arbitrary  criterion  for  rejecting  data 
was  established.  A  significance  level  {oc )  of  o10  was 
chosen  for  the  individual  runs0  This  means  that  we 
accepted  a  10$  chance  of  making  an  error  when  rejecting  a 
result  as  false;  or,  conversely,  only  10$  of  the  time? 
over  a  series  of  many  runs,  would  true  results  fall  outside 
of  our  criterion o   In  the  example  mentioned  above  the 
value  of  0C  for  this  run  would  be  o915s,  meaning  a  91 „ % 
chance  of  being  incorrect  if  the  results  were  rejected 
as  false.  As  an  example  of  a  run  that  was  rejected;  the 
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VAR  x  DF  was  6l065<.  From  the  X2  vs0  degrees  of  freedom 
tables,  we  find  that  there  is  only  a  5%   chance  of  being 
incorrect  in  rejecting  this  result  /14A 
3.   Method  of  Weighting  Results 

Results  of  the  half  life  measurements  were  weighted 
according  to  the  "goodness  of  fit"  of  the  least  squares 
analysis  0 

The  weight  factor,  C*  ,  applied  to  each  result  is: 

1 

where:  Q\   s  (VAR  x  AM"1)* 

VAR  -  weighted  variance  of  fit 
AM     s  the  inverse  of  the  least  squares 
matrix  in  the  analysis „ 
This  G%   value  is  the  value  given  as  the  half  life 
"sigma"  (Appol)  in  the  computer  program  output  for  the 
particular  run  considered., 

Using  this  weight  factor  the  half  life,  ti?  is  found 
by: 

c  I  H 1 


where:  ti-  =  the  half  life  for  the  i   run 
n  s  number  of  runs  usedo 
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The  standard  deviation  of  the  half  life,  (?Mti)8  is 


calculated  from: 


CMti)  * 


n 

z  0i 

l"    1 


i 
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V.   Results 

In  the  analysis  of  the  data  the  mode  of  decay  (i0eos  two 
component  or  one  component)  used  in  the  program  depended  on  the 
half  life  being  measured .   In  the  case  of  Hf^^m  a  two 
component  decay  was  assumed „  For  Hf^7°m  and  Yb-*-77m  Single 
component  decays  were  usedo  This  required ,  in  the  case  of 
Hfl79m?  that  the  shorter  half  life  be  allowed  to  decay  to  a 
small  value  (less  than  1%   of  total)  before  measurements  were 
taken. 

The  tabulated  results  of  the  individual  acceptable  runs 
for  the  measurements  of  the  half -lives  of  the  isotopes  under 
consideration  are  as  follows : 

Hf177  (n,  &M  Hf178m  (,326  Mev  transition  measured) 

Half-life  Sigma 

3.774  sec,  ,177  seco 

3c 996  o209 

3.952  O057 

3.917  o060 

Hf17^  (n, ^ )  Hf179m  (.217  Mev  transition  measured) 

Half-life  Sigma 

13.435  sec,  .092  sec0 

13.519  .073 

13.509  .163 

13.574  o099 

13.373  o090 

Yb176  (n,  tf )  Yb177m  (.212  Mev  transition  measured) 

Half-life  Sigma 

6.564  sec,  o092  seco 

6.470  ,046 

6.466  0O46 

60 542  o091 
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The  weighted  averages  resulting  from  these  measurements  ares 

Isotope  Half-life 

HfWm  (3o92  t   005)seco 

Hf}79m  (1S049  ±   .04)sec0 

Yb^ym  (6o49  t   O03)seco 

The  value  of  1S.49  sec  for  the  half -life  of  Hfl?9m  agrees 
with  the  value  of  IS e  6  t   *#  reported  by  Alexander  and 
Brinckmanne 

The  value  of  3o92  sec0  for  the  half -life  of  Hf^m 
falls  almost  exactly  between  the  measurements  of  3o5  sec„ 
reported  by  Campbell  and  Fettweis,  and  the  4°3  sec  value 
reported  by  Alexander  and  Brinckmann0 

The  value  of  6„49  sec0  for  the  half-life  of  Yb177m 
is  in  good  agreement  with  the  value  of  6o5  secs  also  reporter 
by  Campbell  and  Fettweis0 

Assuming  that  there  are  no  appreciable  systematic 
errors ,  it  would  appear  that  the  present  system  of 
recording  and  analysing  decay  data  is  capable  of  good 
precision .   Thus  we  are  led  to  prefer  the  single  channel 
device  for  the  following  reasons? 

1)  The  system  dead  time  is  considerably  less  than  that 
encountered  with  multichannel  analysers  (lo4  micro-sec  vsc 
20-50  micro-sec  for  many  multi-channel  devices )0   Smaller 
dead  time  allows  working  with  higher  activities 9   hence 
better  precision  in  counting  data* 
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2)  The  dead  time  factor  and  overall  timing  error  can 
be  accurately  determined  for  the  single  channel  instrument. 
Variations  which  occur  in  "live"  (on)  times  and  switching 
(off)  times  between  channels  make  evaluation  of  actual 
counting  times  and  dead  times  a  difficult  task  with  multi- 
channel devices . 

It  is  felt  that  the  present  system ,  with  the  FRANTIC 
computer  program,  can  be  effectively  used  to  measure  half- 
lives  on  the  order  of  0o5  seconds 0 
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Appendix  II 
Isotope  Energy  Levels  and  Decay  Schemes 
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APPENDIX  III 

Isotopic  Analysis  of  Samples 

lo  Analysis  of  Hafnium  Sample 
(as  reported  by  Oak  Ridge  Analysis) 


Element  Hafnium 
Isotope    177 
Series     G.  P. 
Sample     827  (a) 


Isotopic  Analysis 


Atomic 
percent 
trace  on 

ly 

Precision 

o99 

* 

o07 

59,08 

♦ 

o52 

31.73 

+ 

o49 

4,53 

+ 

,17 

3o68 

+ 

oil 

Isotope 
174 

I 

177 
178 

197 
180 

The  limits  emoted  above  are  an  expression  of  the 
precision  of  this  measurement  only.   The  error  is  estimated 
at  less  than  1%   from  known  sources  of  systematic  errors 0 

Spectroscopic  Analysis 

Element  Percent 

Cu  Ool 

Fe  0,2 

2,  Analysis  of  Ytterbium  Sample 
(as  reported  by  Oak  Ridge  Analysis  ) 
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Element  Ytterbium 

Isotopic  Analysis 

Isotope      176 
Series       IU 

Isotope 

Atomic 
percent 

Precision 

Sample      1061  (a) 

168 

<   oOl 

170 

<  002 

171 

0  JL  J.. 

172 

o26 

i  °°5 

173 

o30 

i  005 

174 

i.fti 

t  o05 

176 

97o57 

0  X 

The  limits  quoted  above  are  an  expression  of  the  precision 
of  this  measurement  only0   The  error  is  estimated  at  less  than 
1%   from  known  sources  of  systematic  errors „ 

Spectroscopic  Analysis 
Al  <  o05 
Ba  <  .02 
Be  <  .001 
Ca  <  .07 
Co  <  ,05 
Cr  <  .05 
Cu  <  .05 
Fe  <  .02 
K  <  .01 
Li  <  .  01 

Kg  <  .02 

Mn  <  .02 

Mo  <i,05 

Na  <,05 


Ni 

< 

o05 

Gd 

< 

0OO4 

Pb 

< 

.1 

Tb 

< 

o05 

Si 

< 

.05 

Dy 

< 

0  i 

Sn 

< 

.05 

Ho 

< 

o02 

Ti 

< 

o02 

Fr 

< 

0OO4 

V 

< 

.02 

Tm 

< 

o02 

Zr 

< 

.1 

Lu 

< 

oOlg 

Y 

< 

a004 

La   <  oOl 
Ce   <  .2 
Pr   <  .1 


Nd 

< 

ol 

Sm 

< 

ol 

Eu 

< 

35 

.004 

<=  No  spectrum  line 
visible  0  Probably 
absent 3  definitely 
less  than  value  given, 


